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Employing first-principle electronic structure calculations, we study the magnetic and electronic 
properties of the XY (X= Mg, Ca, Sr and Y= N, P, As, Sb) compounds crystallizing in the 
zincblende structure. The Ca and Sr alkaline-earth metal monopnictides are found to be half- 
metallic with a total spin magnetic moment per formula unit of 1.0 [is- In the case of the Mg 
alloys the p-d hybridization effect is much weaker and only MgN is a half-metal. Electron counting 
of the bands explains the Slater-Pauling behavior exhibited by the total spin magnetic moment. 
We also study for these alloys the effect of deformation taking into account both the cases of 
hydrostatic pressure and tetragonalization keeping constant either the in-plane lattice parameters 
or the unit cell volume. Even large degrees of deformation only marginally affect the electronic 
and magnetic properties of these alloys. Finally, we show that this stands also for the rocksalt 
structure. Our results suggest that alkaline-earth metal monopnictides are promising materials for 
magnetoelectronic applications. 
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I. INTRODUCTION 

Half-mctals are among the most studied magnetic 
materials due to their potential application in spin- 
tronic devices These unusual materials present metal- 
lic properties for one spin-direction while they behave 
like semiconductors for the other spin-direction eventu- 
ally leading to 100% spin-polarization at the Fermi level4 
This peculiar property is exhibited by several transition- 
metal alloys but lately half-metallic ferromagnetism has 
been predicted for several alloys which do not include 
transition-metal atoms. Such alloys are widely known 
in literature with various names like d magnets, p- 
ferromagnets or sp-electron ferromagnets^ These com- 
pounds have the advantage of being energy-efficient for 
applications since they create weak external magnetic 
fields and thus lead to minimal energy losses. There are 
several ways to create sp-electron ferromagnets and an 
extensive review is given in Ref. 0. 

One of the most promising routes to half-metallic sp- 
electron ferromagnets is the growth of I/IFIV/V nanos- 
tructures in metastablc lattice structures similar to the 
case of transition-metal pnictides and chalcogenides in 
the metastablc zincblende structured Several studies to 
this research direction have appeared following the pio- 
neering papers published by Geshi et al 7 - and Kusakabe 
et al^ who have shown using first-principles calculations 
that CaP, CaAs and CaSb alloys present half-metallic 
ferromagnetism when grown in the zincblende structure. 
Moreover they have shown that the gap appears due to 
a p-d hybridization mechanism which is similar as we 



will discuss in the next paragraph for all half-metallic 
ferromagnetic I/IFIV/V alloys in all three zincblende, 
wurtzite and rocksalt mctastable structures. The spin 
magnetic moment in these alloys follows a Slater-Pauling 
behavior with the total spin magnetic per formula unit 
in hb being 8 minus the number of valence electron in 
the unit cell: M t = 8 — Z t . The eight stems from the 
fact that in the majority-spin band we have exactly 4 oc- 
cupied states (one s and three p states) and the rest of 
the electrons occupy exclusively the minority-spin states. 
The same reasoning stands even for Hcuslcr compounds 
like GeKCa and SnKCa which contain no transition 
metals and have been predicted to exhibit half-metallic 
ferromagnetism^ Evidence of the growth of such nanosc- 
tructures has been provided by Liu et al who have re- 
ported successful self-assembly growth of ultrathin CaN 
in the rocksalt structure on top of Cu(001)J£ 



II. OVERVIEW OF RECENT LITERATURE 

Following the Rcfs. and[H mentioned in the previous 
paragraph a lot of theoretical studies on sp-electron ferro- 
magnets have appeared based on first principles calcula- 
tions and we will give a short overview of them. A general 
trend is that the rocksalt (RS) structure is more stable 
than the zincblende (ZB) or wurtzite (WZ) structures, in 
all cases the Slater-Pauling (SP) rule for the total spin 
magnetic moment is obeyed and the spin magnetic mo- 
ment is almost entirely concentrated to the anionic sites 
and the interstitial region while the cations show a spin 



moment one or two orders of magnitude smaller (note 
that the occupied states are almost exclusively of an- 
ionic s- and p-charactcr). Also in most studied cases 
half-metallicity is conserved for large variations of the 
lattice constant assuming hydrostatic pressure. 

We will start our overview from the alkali-metal based 
alloys. In 2006 Sicberer and collaborators predicted half- 
metallicity for LiN and Na(K)P(As) in the zincblende 
structure with a total spin magnetic moment of 2 [Ib 
per formula unit (f.u.)^ and in 2008 Zhang confirmed 
the half-metallicity also for the alkali-metal carbides 
(Li,K,Na)C with a total spin moment of 3 /is in ac- 
cordance to the SP ruled After. Zberecki and collab- 
orators have demonstrated that the rocksalt is the most 
stable structure for alkali-metal nitrides with respect to 
the zincblende and wurtzite structures^ These results 
were independently confirmed by Yand and by Gao and 
collaborators who have also shown that among the alkali 
metal phosphides and arsenides, only KP and KAs are 
half- metallic fcrromagnets. ls Finally, Gao et al have also 
studied the alkali-metal sulfides (Li,Na,K)S which exhib- 
ited a total spin magnetic moment of 1 [Ib as expected by 
the SP rule for both the ZB and the energetically-favored 
RS structures d 

Gao and collaborators have studied a series of alkaline- 
earth metal compounds with the IVth column elements 
and found that all CaSi, CaGe, CaC, SrC and BaC in 
the ZB structure are half-met alsi 17 ' 18 Gao and Yao have, 
then, studied SrC and BaC in more detail and have shown 
that the RS structure in the most stable against the other 
metastable structures and found that both RS-SrC and 
RS-BaC are half-metals with the Curie temperature of 
SrC exceeding the room temperature (485 K)d Dong 
and Zhao have shown that the alkaline-earth metal car- 
bides retain their half-metallicity up to a hydrostatic 
pressure of 32GPa (the volume is reduced by 30%),— 
Zhang et al have studied ScC in the WZ structured and 
finally Gao have shown that the (110) surface of ZB-CaC 
keep their half- metallic character contrary to the (001) 
surfaces d Also Verma and collaborators have studied 
the alkaline-earth metal silicides and found that only the 
ZB-structure is half-metallic contrary to the more-stable 
RS ferromagnetic structured 

Most of the attention has been focused on the alkaline- 
earth metal (Ilnd column) compounds with the Vth- 
column elements. As we have already mentioned half- 
metallicity has been predicted for the CaP, CaAs and 
CaSb ZB- alloys Sicberer et al studied all possible II- 
V combinations in the ZB and WZ structures and found 
that all alloys containing Ca, Sr and Ba are half-metallic 
while only MgN was half-metallic between the Mg-bascd 
compounds d It was also shown in Ref. [ll|that the ferro- 
magnetic state is energetically preferable to both the non- 
magnetic and the antiferromagnetic configurations. Yao 
et al have shown that the ZB-alloys containing Bi from 
the Vth column are also half metalsd but later on Li and 
Yu demonstrated that the inclusion of spin-orbit coupling 
destroys half-metallicity in these Bi-alloysd Volnianska 
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FIG. 1: (Color online) Calculated equilibrium lattice con- 
stant in Angstr0m of several alkaline-earth pnictides assum- 
ing the zincblende structure. Total energy calculations have 
been performed using both the local-spin density (LSDA) and 
generalized gradient (GGA) approximations to the exchange- 
correlation potential in conjunction to the full-potential 
nonorthogonal local-orbital minimum-basis band structure 
method (FPLO). 



and Boguslawski, as well as Geshi and collaborators have 
studied the alkaline-earth metal nitrides and have shown 
that the RS is the most stable structure with formation 
energies of about -11 eV per unit cel L 26 ' 27 Gao et al have 
shown that among the RS alloys containing Ca, Sr or Ba 
as a cation and N, P or As as an anion only the nitrides 
are stable half-metallic ferromagnets with a total spin 
magnetic moment of 1 \ib and cohesive energies about -9 
eV per formula unit.— Droghetti and collaborators have 
shown that RS-MgN is in verge of the half-metallicity 
and suggested that MgN inclusion upon the N-doping of 
MgO should lead to a material suitable for magnetic tun- 
nel junctionsd The most recent studies on nitrides con- 
cern the Curie temperature which was found to be well 
above the room temperature for CaN and SrN in both 
the ZB and RS structures^! and the RS-CaN/ZB-InN 
and RS-SrN/ZB-GaP (111) interfaces which were found 
to retain half-metallicity only when the interface is made 
up from Ca-N or N-In atoms in the first case and N-Ga 
in the second cased 



III. MOTIVATION AND COMPUTATION 
METHOD 

From the discussion in Section |TT] we can conclude 
that alkaline-earth metal pnictides are very promising for 
spintronic applications since (i) have a small spin mag- 
netic moment per formula unit (1 hb) and thus create 
small external magnetic fields, (ii) they present very sta- 
ble half-metallicity upon hydrostatic pressure, (iii) their 
equilibrium lattice constant are close to a lot of semi- 
conductors, (iv) results on the ZB and RS structures 



3 



TABLE I: Equilibrium lattice constants in Angstr0m determined through total energy calculations employing the FPLO method 
within both the LSDA and GGA approximations for all zincblende XY compounds under study. With m x and m Y we present 
the atom-resolved spin magnetic moment in at both the cation (X) and anion (Y) sites at the corresponding equilibrium 
lattice constants, as well as the total spin magnetic moment (m total ) per formula unit. The Mg-based alloys, with the exception 
of MgN, have converged to non-magnetic ground states. 



XY 




(A) 


m 


X (Wb) 






m total (u R ) 


comp ounds 


LSDA 


GGA 


LSDA 


GGA 


LSDA 


GGA 


LSDA GGA 


MgN 


4.74 


4.84 


-0.054 


-0.080 


1.054 


1.080 


1.0 1.0 


MgP 


5.78 


5.89 


Non- 


magnetic 


Non-magnetic 


Non-magnetic 


MgAs 


5.95 


6.09 


Non- 


magnetic 


Non-magnetic 


Non-magnetic 


MgSb 


6.40 


6.56 


Non- 


magnetic 


Non-magnetic 


Non-magnetic 


CaN 


5.31 


5.45 


-0.053 


-0.073 


1.053 


1.073 


1.0 1.0 


CaP 


6.39 


6.56 


0.059 


0.025 


0.941 


0.975 


1.0 1.0 


CaAs 


6.55 


6.75 


0.082 


0.044 


0.918 


0.956 


1.0 1.0 


CaSb 


7.03 


7.24 


0.144 


0.102 


0.856 


0.898 


1.0 1.0 


SrN 


5.66 


5.82 


-0.071 


-0.093 


1.071 


1.093 


1.0 1.0 


SrP 


6.74 


6.91 


0.003 


-0.031 


0.997 


1.031 


1.0 1.0 


SrAs 


6.90 


7.10 


0.018 


-0.026 


0.982 


1.026 


1.0 1.0 


SrSb 


7.37 


7.60 


0.065 


0.018 


0.935 


0.982 


1.0 1.0 



suggest high values of the Curie temperature, (v) the 
half-metallic gaps are wide, and (vi) interfaces with semi- 
conductors retain half-metallicity. In this article we will 
study in detail (Section IIV[) the electronic and magnetic 
properties of several pinctides in the zinc-blende struc- 
ture having the chemical formula XY, where the Ilnd 
column element X is Mg, Ca or Sr and the Vth column 
element Y is N, P, As or Sb. This study is followed (Sec- 
tion [Vj by an extensive investigation of the properties 
of these alloys under hydrostatic pressure, tetragonal- 
ization keeping the in-plane lattice parameters constant 
and tctragonalization keeping the volume of the unit cell 
constant. These latter calculations simulate the effect of 
growth on various substrates where one of the three cases 
of deformation can occur. To make our study more com- 
plete we have also conducted calculations for RS CaN 
and SrN under such deformations. 

To perform our calculations we employ the full- 
potential nonorthogonal local-orbital minimum-basis 
band structure scheme (FPLO)^ using a dense 
20x20x20 k-mesh in the Brillouin zone in the recipro- 
cal space to perform integrations. Since we are inter- 
ested in deformations, we have hrst to determine the 
equilibrium lattice constant using total energy calcula- 
tions for all compounds under study in the ZB lattice 
structure. In Fig. [T]we summarize our results using both 
the local-spin density (LSDA)^> and generalized-gradient 
(GGA)^. approximations to the exchange correlation- 
potential. GGA is know to soften the bonds with respect 
to LSDA and to produce equilibrium lattice constants 
more close to experimental data (in fact GGA overesti- 
mates lattice constants by less than 1% while LSDA un- 
derestimates them by more than 1%). Both LSDA and 
GGA curves exhibit exactly the same behavior in Fig. Q] 
with the LSDA lying under the GGA values as expected. 
As we move from Mg to Ca and then to Sr for the same 
Y element the calculated equilibrium lattice constant be- 



comes larger since we move from a lighter to a heavier 
atom. The same trend occurs for one X-element as we 
move from N to P and then to As and Sb. Interestingly 
the phosphides and arsenides have very close lying equi- 
librium lattice parameters. As stated above the GGA is 
expected to give a better estimation of the elastic prop- 
erties with respect to LSDA while both yield similar de- 
scription of the magnetic properties for the same lattice 
constant. In Table fl] we have gathered our results for all 
alloys under study using both LSDA and GGA. Our cal- 
culated GGA equilibrium lattice constants for CaP, CaAs 
and CaSb are in excellent agreement with the GGA re- 
sults by Geshi and collaborators using the full-potential 
linearized augmented plane waves methods (FLAPW)^ 
Notice that for MgP, MgAs and MgSb both approxima- 
tions converged to a non-magnetic solution. For all other 
alloys both LSDA and GGA give a total spin magnetic 
moment of 1.0 fis at the LSDA and GGA equilibrium 
lattice parameters and thus these alloys are expected to 
present half-metallic behavior. The atom-resolved spin 
magnetic moments show a small variation in the second 
decimal digit between GGA and LSDA. Its origin is the 
different equilibrium lattice constant produced by LSDA 
and GGA, and for the same lattice constant both LSDA 
and GGA produce almost identical spin magnetic mo- 
ments. Thus in the rest of the paper all presented re- 
sults have been obtained using the GGA approximation 
at the GGA equilibrium lattice constants. In general our 
results agree well with the previous ab-initio studies on 
these materials discussed in Section ILT1 



IV. EQUILIBRIUM PROPERTIES 

As shown in Table Q] MgP, MgAs and MgSb converge 
to non-magnetic solutions while all other alloys seem to 
be half-metallic ferromagncts as derived from the Slater- 
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FIG. 2: (Color online) Total DOS in the primitive cell for 
all studied compounds for a wide energy spectrum within the 
GGA approximation at the equilibrium lattice constant. The 
zero energy has been chosen to represent the Fermi energy 
(Ef). The positive DOS values correspond to the majority- 
spin (spin-up) electrons and the negative DOS values to the 
minority-spin (spin-down) electrons. The antibonding states 
lying above both the majority-spin and minority-spin gaps 
are mainly of cationic d-character while the states below the 
gaps are mainly of anionic-p character as shown in Fig. [S] 
Note that we do present the anionic s-states located at about 
-10~-12 eV for both spin directions. 



Pauling rule discussed in Section |TT] which will be again 
referred to later in the text. In Fig. [5] we have drawn 
the total density of states (DOS) for all magnetic alloys 
as a function of the energy setting the Fermi level as 
the zero energy. Moreover positive DOS values refer to 
the majority-spin or spin-up electrons and negative DOS 
values to the minority-spin or spin-down electrons. In 
all cases we do not represent the s-valence states which 
are located at about 10^12 eV below the Fermi level. 
Note that for all compounds under study there are very 
large gaps in both majority- and minority-spin electronics 
bands structures. The Fermi level crosses the minority- 
spin states while it falls within the gap for the majority- 
spin states and thus we conclude that these alloys are 
half-metallic ferromagnets where the majority-spin band 
shows semiconducting behavior while the minority-spin 
band shows metallic behavior. As we move from MgN 
to CaN and then to SrN the bands become more narrow 
while as we move along the Vth column (N— >P— >As— 5-Sb) 
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FIG. 3: (Color online) Atom- and spin-resolved density of 
states (DOS) of the four alloys presented also in Fig. f3] We 
focus around the Fermi level and do not present the deep- 
lying s-states (see caption of Fig. [2]) and the high-energy 
antibonding states (see Fig. [2]for a broader energy spectrum). 
Details as in Fig. [2] 



the gap becomes smaller. 

The main weight of the states lie above the Fermi level. 
This is easy to understand if we examine the origin of 
the gap discussed in Refs. and H. The X-elements 
from the Ilnd column (Mg,Ca,Sr) provide two valence 
electrons (occupying e.g. in Ca the 4s states in the free 
atom) while the anions (N,P,As,Sb) provide 5 valence 
electrons (e.g. in free atom of As the atomic configura- 
tion is 4s 2 4p 3 ). In total there are 7 valence electrons 
per unit cell. The first two occupy the s-valence states 
created by the anions which lie deep in energy as men- 
tioned above. The p-states of anions (N,P,As,Sb) hy- 
bridize strongly with the unoccupied triple-degenerated 
t2g d-states of the cation, which transform following the 
same symmetry operations in case of lattice with tctrahc- 
dral symmetry like the zincblcnde one, and form bonding 
and antibonding hybrids which are separated by large 
energy gaps. The bonding hybrids contain mostly p- 
admixturc while the antibonding hybrids are mainly of d- 
character. The remaining 5 valence electrons occupy the 
bonding hybrids which are mainly of anionic p-character 
in such a way that all three majority-spin p-states are 
occupied while in the minority-spin band the Fermi level 
crosses the bands so that only the two out of three p- 
states are occupied. This gives in total a spin magnetic 
moment per formula unit of exactly 1 /is. In the case 
of Ca and Sr the unoccupied 3d and 4d states, respec- 
tively, are relatively close in energy to the valence p-states 
of the anion and thus hybridization is strong enough so 
that all resulting compounds are magnetic. The case of 
Mg should be discussed separately since for light Mg the 




FIG. 4: Spin-resolved band structure of MgN, CaN, SrN and CaP compounds within the GGA approximation along the high- 
symmetry lines in the Brillouin zone. The zero energy has been chosen to represent the Fermi energy. The solid lines correspond 
to the majority-spin (spin-up) electrons and the dashed lines to the minority-spin (spin-down) electrons. The deep-lying s-band 
is not shown; it is located at 10~12 eV below the Fermi level. 



unoccupied 3d states are very high in energy with re- 
spect to the occupied states compared to Ca for which 
also the first unoccupied d-shcll is the 3d one. As a re- 
sult in MgY alloys the p-d hybridization is very weak 
and the resulting bands are very broad. In MgP, MgAs 
and MgSb the DOS at the Fermi level is weak and these 
alloys prefer to stay non-magnetic and only MgN gains 
energy when it becomes magnetic. Sieberer et al using 
the FLAPW method within GGA have predicted that 
MgP and MgAs are magnetic with very small magnetic 
moments (0.08 and 0.02 fis respectively) instead of non 
magnetic M- The difference should be attributed to the 
fact that, as shown also on Ref. [TJ these compounds 
are at the verge of being magnetic and a small shift of 
the predicted Fermi level or a slightly different shape of 
the predicted DOS by different ab-initio methods could 
eventually lead to the loss of the magnetism,. The lat- 
ter compounds can be also seen as a liming case of the 
Stoner criterion for transition-metal alloys which states 
that magnetism is the ground state only when the DOS 
at the Fermi level in the non-magnetic case is quite high 
and the system gains energy when being magnetic. 

The DOS above the gaps in Fig. [5] comes mainly 
from the unoccupied d-states and exhibits vanishing spin- 



splitting. Also the weight of these states is much larger 
than for the states below the gaps since the latter are 
mainly the p-states (summing to three per spin direction) 
while the former are mainly the d-states (summing to five 
per spin direction). Thus relevant to our discussion on 
the magnetic properties arc mainly the p-d bonding hy- 
brids close to the Fermi level since also the s-states lying 
deep in energy are fully occupied (one electron per spin 
direction). In Fig. [3] we present the atom- resolved DOS 
for the three nitrides under study and CaP. In all cases 
most of the DOS is concentrated on the anions reflect- 
ing the above discussion and the cations carry a much 
smaller weight of the p-d bonding hybrids. The ratio 
between the weight at the cation and the anion site is al- 
most constant among the nitrides but increases slightly 
as we move to the CaP (Ca carries a larger weight of the 
p-d hybrids in CaP than in CaN). The fact that these 
states originate mainly from the p-states of the anions 
is confirmed since the curves for the cations DOS follow 
exactly the shape of the curves for the anions. This be- 
havior of the DOS is easy to interpret since the anionic 
p-states, when expanded around the cationic centers, be- 
come t2 9 -like and thus the cationic d-DOS in the energy 
window presented in Fig. [3] is a weak reflection of the 
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anionic p-DOS. Among the nitrides, CaN and SrN ex- 
hibit similar DOS curves but MgN differs considerably 
since the hybridization effect for this alloy is much less 
intense than for SrN and CaN. Moreover as we substitute 
P for N in CaN the bands become more narrow and the 
exchange-splitting between the majority- and minority- 
spin bands decreases (difference in the energy position 
of the mass centers of the bands). All these features 
can be identified when we look at the band structure for 
these four alloys and in Fig. 2] we have drawn the band- 
structure for both spin-directions in the energy window 
from -5 to +5 eV along the high-symmetry axis in the 
Brillouin zone. Majority-spin (solid lines) and minority- 
spin (dashed lines) bands have exactly the same shape 
for all four compounds and they are simply shifted in en- 
ergy Similarities are clear if we compare CaN and SrN 
where the difference in the energy positions of the bands 
is small. In MgN the bands are much broader while in 
CaP the exchange splitting is smaller and with respect to 
CaN more unoccupied bands appear in the same energy 
window. As a result of the reduced exchange splitting 
also the bonding p-d hybrids are more narrow. Our calcu- 
lated band structures for the Ca-based alloys are in good 
agreement with the ones calculated by Li and Yu using 
a GGA pseudopotential technique^ Finally we should 
note that our band structure confirms that the states be- 
low the gaps are p-d hybrids with mainly p-character, 
since they are triple degenerated at the T point, a char- 
acteristic of either the p or the t2 S states in the lattices 
of tetrahedral symmetry likes the zincblcnde one. 

In the last paragraph we shall also discuss the spin 
magnetic moments presented in Table [I] The cations 
carry very small spin magnetic moments of less than 
0.1 1 fiB and almost all spin moment is carried from the 
anions. The spin magnetic moment reflects directly the 
imbalance in the total number of occupied states. In 
Fig. [5] we have drawn both the DOS and the integrated 
DOS for the bonding p-d hybrids in the case of CaN. 
We do not include in the integrated DOS the deep-lying 
s-statcs which accommodate one electron per spin direc- 
tion and thus carry no net spin moment. If we integrate 
the total DOS we have three states per spin direction 
as expected by the band-structure discussed just above. 
All the majority-spin states are occupied while the in- 
tegrated DOS line cross the Fermi level at the value of 
two. Thus in the minority-spin band only the two out 
of three states are occupied and the total spin magnetic 
moment is 1 [1b since we have an imbalance in the elec- 
tronic charge counting of one. This picture agrees with 
the Slater-Pauling behavior for the sp-electron ferromag- 
nets which we have discussed in Section HI If we now look 
at the Ca resolved DOS in total we have slightly less than 
0.5 electrons in the majority-spin band and about 0.65 
in the minority-spin band. All majority states are oc- 
cupied while in the minority-spin band only the 0.5 out 
of 0.65 arc occupied. This small imbalance in the elec- 
tron counting leads to a small negative spin magnetic 
moment at the Ca site of -0.073 [is- Small variation of 
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FIG. 5: (Color online) Spin-resolved DOS of the p-d bonding 
hybrids (shaded line) and the corresponding integrated DOS 
(thick red line) for CaN within the GGA approximation and 
its projection on the Ca and N atoms. We do not present 
the s-states which are deep in energy and contain one elec- 
tron per spin-direction with their weight almost exclusively 
at the N-sites. Notice that the p-d hybrids can accommodate 
3 electrons per spin direction. All three majority-spin states 
(~2.5 at the N site and ~0.5 at the Ca site) are occupied, 
while only two out of three (~1.5 at the N site and ~0.5 at 
the Ca site) are occupied resulting in a total spin magnetic 
moment of 1 fiB- 



the values of the cationic integrated DOS at the Fermi 
level leads to small variation of the cationic spin mag- 
netic moment around zero. At the N site we have about 
2.5 states in the majority-spin band which arc completely 
occupied while in the minority-spin bands we have about 
2.35 states in total but only ^1.5 are occupied leading 
to a N spin magnetic moment of about 1 /ib- Previous 
calculations led to a similar picture of the atom resolved 
spin-magnetic moments with the cation spin moments 
being one order of magnitude smaller than the anion spin 
magnetic momentSi£i2£ Thus to conclude, we should have 
in all cases a total spin magnetic moment of 1 \ib since 
we have seven valence electrons, all compounds are half- 
metals and the atom resolved spin moments fluctuate 
around or 1 [lb for cations and anions, respectively, 
depending on the p-d hybridization intensity. 



V. DEFORMATIONS 

To simulate deformation of the lattices with respect to 
equilibrium we took into account three cases: (i) hydro- 
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FIG. 6: (Color online) Total DOS of CaN (left panels) and CaP (middle panels) compounds in the zincblende structure and 
CaN (right panels) in the rocksalt structure within GGA at the equilibrium lattice constants compared to : (upper panel) 10% 
uniform expansion or contraction of the lattice as in the case of hydrostatic pressure, (middle panel) tetragonalization keeping 
the in-plane lattice parameters equal to the equilibrium ones and expanding or contracting the out-of-plane lattice parameter 
by 10%, and (lower panel) tetragonalization expanding or contracting the in-plane lattice parameters by 10% and in the same 
time changing accordingly the out-of-plane lattice parameter to keep the unit cell volume constant to the equilibrium one. 



TABLE II: Calculated spin magnetic moments (in fis) for CaN and CaP in the zincblende structure and CaN in the rocksalt 
structure using the FPLO method within the GGA approximation. The first line corresponds to the ideal equilibrium lattice 
constants and we have considered three cases of deformation: (i) hydrostatic pressure (Hydro) which is a uniform compression 
and expansion of the cubic lattice and the percentage in the second column shows the change in the lattice constant ("-" means 
compression and "+" means expansion), (ii) tetragonalization of the lattice keeping the in-plane lattice parameters constant 
and equal to the equilibrium values and changing only the lattice parameter along the c-axis (Tetra-I) and the percentage in the 
second column shows this change, and (iii) tetragonalization where we vary the in-plane lattice parameters by the percentage 
shown in the second column and in the same time we vary also the lattice parameter along the c-axis so that the unit-cell 
volume is kept equal to the equilibrium one (Tetra-II). In all cases the total spin-magnetic is kept equal to 1.0 [ib (even in 
the case of ±10% hydrostatic pressure where the unit-cell volume changes by almost 30% ) and thus the half-metallicity is 
preserved. 

CaN-ZincBlende (5.45 A) CaP-ZincBlende (6.56 A) CaN-RockSalt (5.02A) 



Case 


Ca 


N 


Total 


Ca 


P 


Total 


Ca 


N 


Total 


Ideal 


-0.073 


1.073 


1.0 


0.025 


0.975 


1.0 


-0.065 


1.065 


1.0 


Hydro 


-1% 


-0.072 


1.072 


1.0 


0.025 


0.975 


1.0 


-0.061 


1.061 


1.0 




-5% 


-0.071 


1.071 


1.0 


0.025 


0.975 


1.0 


-0.059 


1.057 


0.999 




-10% 


-0.068 


1.068 


1.0 


0.026 


0.974 


1.0 


-0.048 


1.017 


0.969 




+1% 


-0.074 


1.074 


1.0 


0.025 


0.975 


1.0 


-0.067 


1.067 


1.0 




+5% 


-0.075 


1.075 


1.0 


0.026 


0.974 


1.0 


-0.070 


1.070 


1.0 




+10% 


-0.080 


1.080 


1.0 


0.026 


0.974 


1.0 


-0.079 


1.079 


1.0 


Tetra-I 


-1% 


-0.072 


1.072 


1.0 


0.025 


0.975 


1.0 


-0.063 


1.063 


1.0 




-5% 


-0.072 


1.072 


1.0 


0.025 


0.975 


1.0 


-0.062 


1.062 


1.0 




-10% 


-0.071 


1.071 


1.0 


0.026 


0.974 


1.0 


-0.058 


1.055 


0.997 




+1% 


-0.073 


1.073 


1.0 


0.025 


0.975 


1.0 


-0.065 


1.065 


1.0 




+5% 


-0.073 


1.073 


1.0 


0.026 


0.974 


1.0 


-0.066 


1.066 


1.0 




+10% 


-0.074 


1.074 


1.0 


0.027 


0.973 


1.0 


-0.068 


1.068 


1.0 


Tetra-II 


-1% 


-0.073 


1.073 


1.0 


0.026 


0.974 


1.0 


-0.064 


1.064 


1.0 




-5% 


-0.072 


1.072 


1.0 


0.027 


0.973 


1.0 


-0.064 


1.064 


1.0 




-10% 


-0.071 


1.071 


1.0 


0.027 


0.973 


1.0 


-0.063 


1.063 


1.0 




+1% 


-0.073 


1.073 


1.0 


0.026 


0.974 


1.0 


-0.064 


1.064 


1.0 




+5% 


-0.072 


1.072 


1.0 


0.027 


0.973 


1.0 


-0.064 


1.064 


1.0 




+10% 


-0.071 


1.071 


1.0 


0.029 


0.971 


1.0 


-0.062 


1.062 


1.0 



8 



static pressure (Hydro) which is a uniform compression 
and expansion of the cubic lattice, (ii) tetragonalization 
of the lattice keeping the in-plane lattice parameters con- 
stants and equal to the equilibrium values and changing 
only the lattice parameter along the c-axis (Tetra-I) , and 
(iii) tetragonalization where we vary both the in-plane 
and out-of-plane parameters keeping the unit cell vol- 
ume constant to the equilibrium (Tetra-II). In Table [TT1 
we have gathered our results concerning the spin mag- 
netic moments for CaN and CaP adopting the ZB-lattice. 
We have also performed similar calculations for SrN and 
SrP but the behavior of the Sr compounds is identical to 
the Ca ones and thus we do not present them. Moreover 
we present results also for CaN in the rocksalt structure 
which has been also widely studied in literature. Prelim- 
inary results on SrN in both ZB and RS structure taking 
into account only the case of tetragonalization keeping 
the unit cell volume constant have been published in Rcf . 
l3ll . The latter type of distortion is expected to occur 
more often when these compounds are grown on top of 
various substrates, but the other two cases under study 
arc also susceptible to occur and we include them aim- 
ing to give a better description of the properties of these 
alloys under various distortions. In the second column 
of the table we present for each one of the three cases 
the percentage of change; "-" means compression by that 
percentage and "+" means expansion. Moreover in the 
case of the hydrostatic pressure the percentage refers to 
the lattice constant since the expansion or compression is 
uniform, in the case of Tetra-I where we keep constant the 
in-plane lattice parameters the percentage refers to the 
out-of-plane lattice parameter, and finally in case Tetra- 
II where we keep the unit-cell volume constant the per- 
centage refers to the in-plane lattice constants. For each 
case we took six values of the percentage ±1 %, ±5 % 
and ±10 %. The larger changes in the unit cell volume 
occur for the Hydrostatic case where ±10 % change of 
the lattice constant means compression or expansion of 
the unit cell volume by about 27 %. Amazingly even for 
these extreme cases CaN and CaP in the ZB lattice re- 
main half-metallic with a total spin magnetic moment of 
1 [Ib- CaN in the RS-lattice shows similar behavior to 
the ZB one and only in the case of 10% compression it 
is at the verge of half-mctallicity and Fermi level is lo- 
cated at the edge of the occupied majority-spin states; 
this behavior occurs because the half-metallic gap (en- 
ergy difference between the upper edge of the occupied 
majority-spin states and the Fermi level) is 0.1 eV much 
smaller than the compounds in the ZB lattice. Defor- 
mations in CaN lead to small changes of the absolute 
values of the Ca and N atomic spin magnetic moments 
of less than 0.07 [is in such a way that they cancel each 
other keeping the total spin moment constant. When 
we replace P for N in CaN the fluctuations of the atomic 
spin magnetic moments of the constituent atoms are even 
smaller. 

In Fig. |5] we present for the alloys in Table [TT| the to- 
tal DOS taking into account expansion or compression 



by 10% which are the limiting cases studied and com- 
pare them with the DOS for the ideal cubic lattices at 
the equilibrium lattice constants. We focus on the en- 
ergy window including the p-d bonding hybrids which 
are the most relevant for magnetism. In the case of the 
±1% and ±5% cases the DOS for the distorted lattices 
coincide with the DOS for the equilibrium lattice and the 
curves fall almost one on top of the other; the only no- 
ticeable difference concerns the hydrostatic pressure case 
where compression leads to slightly wider bands and thus 
smaller pick intensities. Differences in the DOS are only 
visible in the case of the ±10% deformations presented in 
the figure. The most important changes occur in the case 
of hydrostatic pressure where compression leads to wider 
bands and expansion to narrower bands, this behavior 
being typical for p-orbitals. In the case of tetragonaliza- 
tion the ZB-CaN shows very small changes in the DOS 
shape while for ZB-CaP we remark a small redistribution 
of the weight of the empty minority-spin states just above 
the Fermi level. For CaN in the RS-lattice the changes in 
the DOS are even more marginal than for the ZB-lattice. 
The large half-metallic gaps especially in the case of the 
ZB lattice ensure that any small change in the character- 
istics of the DOS does not affect the half-metallicity and 
the Fermi level still crosses exclusively the minority-spin 
bands. The same image is also valid for the SrN and SrP 
alloys. Thus we can safely conclude that deformation 
of the lattice only marginally affects the electronic prop- 
erties as was the case also for the magnetic properties 
discussed just above. 



VI. CONCLUSIONS 

We have studied the properties of the XY (X= Mg, 
Ca, Sr and Y= N, P, As, Sb) pnictides crystallizing in 
the zincblende structure using a full-potential electronic 
structure method in conjunction to the generalized gra- 
dient approximation. We have shown that the Ca- and 
Sr-based alloys should be classified as sp-electron fer- 
romagnets (also known as d°-ferromagnets). Moreover 
they are half-metals since the Fermi level crosses only 
minority-spin states and the total spin magnetic moment 
per formula unit was found to be 1.0 [Ib- Contrary to 
the Ca and Sr alloys, in the case of the Mg alloys the 
p-d hybridization effect is much weaker and only MgN is 
a half-metallic ferromagnet while MgP, MgAs and MgSb 
prefer energetically to be non-magnetic. We have also 
studied for these alloys the effect of deformation taking 
into account both the cases of hydrostatic pressure and 
tetragonalization keeping constant either the in-plane lat- 
tice parameters or the unit cell volume. Even large de- 
grees of deformation only marginally affected the elec- 
tronic and magnetic properties of these alloys which kept 
their half-metallic ferromagnetic character. Similar be- 
havior is found also when they crystallize in the rocksalt 
instead of the zincblende structure. Our results suggest 
that the alkaline-earth metal monopnictides arc promis- 
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ing materials for spintronic applications. 



Electronic address: kozdogan@yildiz.edu.tr 

Electronic address: galanakis@upatras.gr 

I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 

323 (2004). 

C. Felser, G. H. Fecher, and B. Balke, Angew. Chem. Int. 
Ed. 46, 668 (2007). 

H. Zabel, Materials Today 9, 42 (2006). 
M. I. Katsnelson, V. Yu. Irkhin, L. Chioncel, A. I. Lichten- 
stein, and R. A. de Groot, Rev. Mod. Phys. 80, 315 (2008). 
O. Volnianska and P. Boguslawski, J. Phys.:Condens. Mat- 
ter 22, 073202 (2010). 

Ph. Mavropoulos and I. Galanakis, J. Phys. Condens. Mat- 
ter 19, 315221 (2007). 

M. Geshi, K. Kusakabe, H. Tsukamoto, and N. Suzuki, 



2004 Preprint arXiv:cond-mat/0402641 (2004) 



K. Kusakabe, M. Geshi, H. Tsukamoto, and N. Suzuki, J. 
Phys. .-Condens. Matter 16, S5639 (2004). 
J. Chen, G. Y. Gao, K. L. Yao, and M. H. Song, J. All. 
Cmps. 509, 10172 (2011). 

X. Liu, B. Lu, T. Iimori, K. Nakatsuji, and F. Komori, 
Surf. Set. 602, 1844 (2008). 

M. Sieberer, J. Redinger, S. Khmelevskyi, and P. Mohn, 
Phys. Rev. B 73, 024404 (2006). 

C.-W. Zhang, J. Phys. D.: Appl. Phys. 41, 085006 (2008). 
K. Zberecki, L. Adamowicz, and M. Wierzbicki, Phys. St. 
Solidi (b) 246, 2270 (2009). 
E. Yan, Physica B 407, 879 (2012). 

G. Y. Gao, K. L. Yao, Z. L. Liu, Y. Min, J. Zhang, S. 
W. Fan, and D. H. Zhang, J. Phys.: Condens. Matter 21, 
275502 (2009). 

G. Y. Gao, K. L. Yao, M. H. Song, and Z. L. Liu, J. Magn. 
Magn. Mater. 323, 2652 (2011). 

G. Y. Gao, K. L. Yao, Z. L. Liu, J. L. Jiang, L. H. Yu, and 
Y. L. Shi, J. Phys.: Condens. Matter 19, 315222 (2007). 



1S G. Y. Gao, K. L. Yao, E. §a§ioglu, L. M. Sandratskii, Z. 
L. Liu, and J. L. Jiang, Phys. Rev. B 75, 174442 (2007). 

19 G. Y. Gao and K. L. Yao, Appl. Phys. Lett. 91, 082512 
(2007). 

20 S. Dong and H. Zhao, Appl. Phys. Lett. 98, 182501 (2011). 

21 C.-W. Zhang, S.-S Yan, and H. Li, Phys. St. Solidi (b) 245, 
201 (2008). 

22 G. Y. Gao and K.-L. Yao, J. Appl. Phys. 106, 053703 
(2009). 

23 U. P. Verma, Mohini, P. S. Bisht, and P. Jensen, Semicond. 
Set. Technol. 25, 105002 (2010). 

24 K. L. Yao, J. L. Jiang, Z. L. Liu, and G. Y. Gao, Phys. 
Lett. A 359, 326 (2006). 

25 Y. Li and J. Yu, Phys. Rev. B 78, 165203 (2008). 

26 O. Volnianska and P. Boguslawski, Phys. Rev. B 75, 
224418 (2007). 

27 M. Geshi, K. Kusakabe, H. Nagara, and N. Suzuki, Phys. 
Rev. B 76, 054433 (2007). 

28 G. Y. Gao, K. L. Yao, Z. L. Liu, J. Zhang, Y. Min, and S. 
W. Fan, Phys. Lett. A 372, 1512 (2008). 

29 A. Droghetti, N. Baadji, and S. Sanvito, Phys. Rev. B 80, 
235310 (2009). 

30 A. Laref, E. §asioglu, and I. Galanakis, J. Phys.: Condens. 
Matter 23, 296001 (2011). 

31 K. Ozdogan, E. §a§ioglu, and I. Galanakis, submitted; 
preprint larXiv:1203.1106l (2012L 

32 G. Y. Gao, K. L. Yao, and N. Li, J. Phys.: Condens. Matter 
23, 075501 (2011). 

33 K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999). 

34 J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 

35 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
78, 1396 (1997). 



